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We constrain energy spectra of supernova neutrinos through the avoidance of an overproduction 
of the n B abundance during Galactic chemical evolution. In supernova nucleosynthesis calculations 
with a parametrized neutrino spectrum as a function of temperature of ^ MlT and and total 
neutrino energy, we find a strong neutrino temperature dependence of the 11 B yield. When the 
yield is combined with observed abundances, the acceptable range of the Vp,, T and P P>T temperature 
is found to be 4.8 to 6.6 MeV. Nonzero neutrino chemical potentials would reduce this temperature 
range by about 10% for a degeneracy parameter r\ v — fi u /kT v smaller than 3. 

PACS numbers: 26.30.+k, 97.60.Bw, 25.30.Pt 



(N 
> 
CO 

O 
in 
o 
in 
o 

On; 

6 

03 



The light elements (Li-Be-B) are continuously pro- 
duced by supernovae (SNe) QHMOl, as well as interac- 
tions of Galactic cosmic rays (GCRs) with the interstel- 
lar medium (ISM), nucleosynthesis in asymptotic giant 
branch (AGB) stars, and novae during Galactic chemical 
evolution (GCE, i.e., the evolution in chemical composi- 
tion of stars and interstellar gas during Galactic history) 
0,H|. In the case of boron, cosmic ray induced spallation 
in the ISM and supernova ejecta dominate the produc- 
tion; n B is contributed through both channels, while 10 B 
production is probably exclusively due to GCRs. The 
contribution from supernovae to the production of n B 
can be calibrated with the isotopic ratio N( 11 B)/N( 10 B), 
measured with great precision in primitive meteorites. 

The SN ^-process plays an important role for 11 B and 
7 Li production nj. The interaction of neutrinos, emitted 
in copious amounts during core collapse and the subse- 
quent cooling phase of proto- neutron stars, with matter 
in the ejecta of SNe, contributes uniquely to GCE. Re- 
cent studies based on the theoretical yields derived by 
Woosley and Weaver (WW95) @ suggest that the SN 
contribution to the n B abundance is significantly larger 
than that required to explain the boron evolution in the 
Galactic disk and the meteoritic n B/ 10 B ratio HHEl. 

To match the abundance of 11 B established during 
GCE, we previously assumed neutrino energy spectra to 
resemble Fermi-Dirac (FD) distributions with zero chem- 
ical potentials \x v = Q, 12, an d fixed neutrino tem- 
peratures of 6.0, 3.2, and 5.0 MeV for ^ MjT (v^t), ^e, an d 
j> , respectively 0. The v^^ T temperature of 6.0 MeV is 
significantly smaller than the 8.0 MeV used in the other 
previous studies of the ^-process 0, 0, ■ This reduc- 
tion is derived from an investigation of the dependence 
of the n B yield on the total neutrino energy E v and the 
decay time scale t„ of the neutrino flux. The yield is 



roughly proportional to E v and rather insensitive to t v . 
The temperature dependence was not investigated very 
well. 

Studies of supernova explosions with detailed neutrino 
transport (e.g., [il Il0| and references therein) have indi- 
cated that emerging neutrino spectra do not closely fol- 
low FD distributions with \x v = 0. Since the high-energy 
tail of the energy distribution is predominantly impor- 
tant for the z/-process (e.g. pj), the use of FD distribu- 
tion with /i„ = may be justified as an approximation as 
long as the spectrum above e v s» 10 MeV is a good match 
to the shapes obtained in detailed transport simulations 
dHUHH^ ■ However, if the 11 B yield depends strongly on 
the neutrino temperatures, which have not yet been clar- 
ified theoretically, the nonzero chemical potentials would 
change the resultant n B abundance in a different matter 
from what follows from FD distributions with /z„ = 0. 
The purpose of this Letter is to investigate the neutrino 
temperature dependence of the SN ^-process in detail, 
and to find out how robustly lower neutrino tempera- 
tures may provide the means to avoid overproduction of 
the n B abundance in GCE and meteoritic 11 B/ 10 B ratio. 

The adopted model for SN neutrinos is guided by nu- 
merical simulations from the literature, with a few addi- 
tional simplifying assumptions. The neutrino luminosity 
is assumed to be uniformly partitioned among the neu- 
trino flavors, and is assumed to decrease exponentially in 
time with a time scale t v = 3 s Q. The latter assump- 
tion is not critical, because the ejected masses of n B and 
7 Li are insensitive to tv . We initially assume that the 
spectra indeed obey the FD form with \x v = 0. 

Only the total neutrino energy E v and the tempera- 
tures T Vfi T are free parameters. The allowed range of the 
total neutrino energy E v is 



1 x 10 53 ergs <E v <6x 10 s3 ergs 



53 , 



(1) 
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which includes the reduced range 2.4 x 10 53 ergs < E v < 
3.5 x 10 53 ergs, corresponding to the estimated range 
in gravitational binding energy of a neutron star with 
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FIG. 1: Mass fraction distribution of the light elements in 
the 16.2M model with T Vfl T = 6 MeV. The horizontal axis 
is the interior mass in units of the solar mass. 



mass ~ IAMq |lj|. The considered range of the neutrino 
temperature T v is 



4.0 MeV < T v ^ < 9.0 MeV. 



(2) 



Temperatures of the v c and D c , T„ o and T Pc , are less im- 
portant for the ^-process of the light elements, and we 
set their values to 3.2 and 5.0 MeV, respectively Q. 

The SN model used in this work is identical to that 
described by 0. We use progenitor model 14E1, with a 
mass at explosion of 16. 2Mq ^4j, corresponding to SN 
1987A. The propagation of a shock wave during the SN 
explosion is followed with a spherically symmetric La- 
grangian PPM code 0, 0] . The explosion energy and 
the mass cut are set to 1 x 10 51 ergs and 1.61Mq, re- 
spectively. Then, we calculate explosive nucleosynthesis 
by postprocessing as described in • The reaction rates 
of the ^-process are derived by interpolating the loga- 
rithmic values of the cross sections listed in the tables 
of |17| . This setup determines the thermodynamic histo- 
ries of the various mass shells that ultimately constitute 
the supernova ejecta (no fall back), and the ^-process 
yields within the ejecta are then determined through the 
cross sections as soon as the time- and energy-dependent 
neutrino flux is specified. We calculate the yields for a 
parameter grid with 126 points, with steps of 1 x 10 53 
ergs in E v and steps of 0.25 MeV in T v t T . Shown in Fig. 
1 is an example of the produced mass fractions for E v — 
3 x 10 53 ergs and T Vv t = 6 MeV. 

Ratios of the ejected masses of 11 B and 7 Li to those of 
WW95 |2j, defined as the overproduction factor /„, are 
shown in Fig. 2 as contours in the E v vs T v T plane. For 
the WW95 case of T VfiT = 8 MeV and £„ = 3x 10 53 ergs, 
we find ejected 11 B and 7 Li masses of 1.92 x 10 _6 Mq and 
7.37 x 1O _7 M0, which are very close to the yields of 1.85 x 
1O~ 6 M and 6.67 x 1O- 7 M obtained with the S20A 
model of respectively. The n B mass ratio changes 
between 0.038 (lower left corner of Fig 2a) and 2.9 (upper 
right corner) in the assumed ranges of E v [Eq. (Q)] and 
T u [Eq. (J2J]. The mass ratio of 7 Li changes between 
0.039 [lower left corner of Fig 2(b)] and 3.3 (upper right 
corner). Note that dependence on the explosion energy 
and the mass cut is weak. 
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FIG. 2: Contour lines of the overproduction factor /„ for (a) 
11 B and (b) 7 Li in the parameter plane of total neutrino en- 
ergy E v and neutrino temperature T u T (see text for details) . 
The region between the two solid vertical lines indicates the 
energy range relevant for a neutron star of mass ~ 1.4M 
p^|. The point labeled WW95 indicates the specific parame- 
ter values used in Q. In panel (a), the region between the two 
solid contour lines is the range of ejected mass appropriate for 
GCE of 11 B. The shaded region is the part of parameter space 
in which both constraints (GCE yield of n B and neutron star 
binding energy) are simultaneously satisfied. A similar box is 
drawn in (b) for the case of 7 Li. 



We constrain the neutrino temperature T v T by requir- 
ing that overproduction of n B in GCE must be avoided. 
The overproduction factor depends on details of the GCE 
model, and ranges between 0.18 and 0.40 Q. These 
values are obtained by combining the solar 11 B/ 10 B ra- 
tio with a measure of the relative cosmic-ray and su- 
pernova contribution to solar n B. They are shown in 
Fig. 2(a) as two solid lines. If we adopt the 1.4M Q neu- 
tron star energy range mentioned above (2.4 x 10 53 ergs 
< E v < 3.5 x 10 53 ergs 0), we obtain the shaded re- 
gion shown in Fig. 2(a), which implies that the neutrino 
temperature T„ satisfies 



4.8 MeV < T„ < 6.6 MeV. 



(3) 



With the neutrino temperature and total energy con- 
strained by GCE of n B, we can derive a corresponding 
constraint on the 7 Li yield. Figure 2(b) shows the shaded 
region corresponding to the E v -T v r limits of the shaded 
box in Fig. 2(a). This region implies an ejected mass ra- 
tio of 7 Li between 0.19 and 0.43. If n B production is in- 
deed dominated by the contributions from the ^-process, 
the analysis presented above implies a predicted range of 
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yields for 7 Li, which in turn constrains the contribution 
to 7 Li production from AGB stars and novae. 

We note that the smallest value of our allowed range 
for T v T is in fact smaller than the assumed value of 
Tp o = 5.0 MeV. Since the neutrinospheres of v c and P e 
are larger than those of v^ T and T due to charged 
current interactions, the average energy of v v and u e are 
smaller than those of i/„ r and v^, T (e.g., |!j)- Thus, if 
T„ T is indeed smaller than 5.0 MeV, T„ o and T Po should 
also be smaller than 5.0 MeV. 

We also note that neutrino oscillations would raise the 
contribution of electron neutrinos to the n B and 7 Li 
production. If neutrino conversion between v e and i/ MiT 
occurs in the oxygen-rich layer [e.fr large mixing angle 
with 6 13 large (LMA-L) case in 0, 0], the rates of 
charged current reactions such as 4 He(^ G , e~p) 3 He and 
12 C(i/ e ,e~p) 11 C increase, keeping the rates of neutral- 
current reactions unchanged. The yields of 11 B and 7 Li 
would increase by this effect and, thus, lower neutrino 
temperature is favorable for avoiding overproduction of 
n B. Additional constraints derive from r-process nucle- 
osynthesis in neutrino driven winds 

We use a specific stellar mass model of ~ 20M Q to 
demonstrate the sensitivity of the n B and 7 Li yields to 
E v and T v T . This sensitivity can also be applied to 
supernova models with different progenitor masses, be- 
cause the dominant production processes for 11 B and 7 Li 
are the ^-process and a-capture reactions, which are in- 
sensitive to progenitor masses, specifically, the 4 He and 
12 C abundances. In the He-rich layer, 7 Li is produced 
through the reaction sequences 4 He(^, v'p) 3 R(a, 7) 7 Li 
and 4 He(V, i/'n) 3 He(a, 7) 7 Be(n,p) 7 Li. Most of 11 B is pro- 
duced through 7 Li(a, 7 ) n B and 7 Be(a, 7) 11 C(/?+) 11 B, 
or the ^-process 12 C(V, i/p) n B in the oxygen-rich layer 
0,0,0- The dependence on E v and T v of the ejected 
masses of n B and 7 Li solely relates to that of the v- 
process reaction rates. The a-capture rates do not de- 
pend on the neutrino parameters, and the abundances of 
4 He and 12 C are solely determined during the precollapse 
stage. Thus, the ejected masses of 11 B and 7 Li are pro- 
portional to the ^-process reaction rates in accordance 
with the values of E v and T v T . The neutrino spectrum 
might depend on progenitor mass, but the extent of this 
effect has not yet been established. 

Many studies of z/-induced nucleosynthesis assume FD 
distributions with At„ = o[l|l|l[l|||. However, simula- 
tions of neutrino transport in supernova explosions show 
that the energy spectra are better represented by FD dis- 
tributions with nonzero chemical potential [ToL ITTl IT^ | . 
Therefore, we now consider the effect of nonzero chemi- 
cal potentials within a semianalytic model. We assume 
that the energy dependence of the neutrino-matter inter- 
action cross sections can be expressed as a simple power 
law cr(e) = (7Q6 a . The specific case of a = 2 was discussed 
in . Here we extend their discussion to a wider range 
of values for a. We assume that the energy spectra are 
exact FD distributions, specified by values of tempera- 
ture T v and degeneracy parameter r\ v = fx^/kT^, where 
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FIG. 3: The scaling function with nonzero chemical potential 
C a (r) v ) denned by Eq. Q. The particular values of a shown 
are 4, 5, 6, and 7. 



k is the Boltzmann constant. With the following defini- 
tion of a moment function 



1 



x q dx 



2ir 2 (hc) 3 J exp(x — rj u ) + f 



(4) 



the neutrino number density n u (T v , T7 V ) and energy den- 
sity t v (Tv,^) can be expressed as F 2 {rj v )(kT v ) and 
F 3 (r] u )(kT^) i 1 respectively. For a neutrino spectrum 
specified by T v and r\ Vl the average cross section a{T v , rj^) 
is then given by \F a+2 (r) v ) / F 2 (r] v )]a (kT u ) a , and is re- 
lated to the average cross section one would obtain 
from a spectrum with zero chemical potential <7(T„,0) 
as [F a+2 (7 h )F 2 (0)}/[F a+2 (0)F 2 {7 lu )}a(T u ,0). 

The reaction rate for any i/-process reaction under con- 
sideration at a given time t and at a distance r from 
the source is given by X(T v ,T] u ;t) = a((T Vl ri l ,)<\)(T l , 1 r\ v \t), 
where the neutrino number flux is 



(j>{T v ,r) v ;t) 



1 



E„ 



1 



4nr 2 ^p^kT v t v 



exp 



t — r/c 



(5) 



Note that <j)(T v ,r) v ;t) is a function of not only T v but 
also r\ v because the average neutrino energy depends on 
T v and r\ v : the average energy per neutrino is {e v ) = 
(F 3 {r]„)/F 2 (7 h ))kT„ and F 3 (0)/F 2 (0) = 3.1514 for 7 h = 
0. The reaction rate \(T Vl rj v ;t) can then be expressed 
as 



\(T VlVv ;t) = CM\{T Ul 0;t) 
where C a (77„) is the scaling function 

F a +M F 3 (0) 



C a {rj v ) 



F a+2 (0) F 3 {i7 v ) ' 



(6) 



(7) 



and both X(T V , r\ v \ t) and A(T„,0;i) have the same T v 
dependence oc T"^ 1 . 

We now apply this semianalytic model. First, we deter- 
mine the effective power law indices for the total neutral- 
current cross sections on 56 Fe and 58 Fe from the calcu- 
lations presented in [2(j. We find indices of 3.7 and 3.8, 
respectively. This implies a(T v ,3)/a(T v ,0) =1.72(1.73) 
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for a=3.7(3.8), consistent with the values reported in 
[20|. We also evaluate the power law indices a of the 
cross sections of 4 He(V, i/p) 3 H and 12 C{v, z/p) n B by fit- 
ting the cross sections in ^7|- For 4 He(V, v'p) 3 K and 
12 C(V, i/p) n B, we find a — 6.7 and 5.9, respectively. 
These indices are much larger than the values obtained 
for reactions with the larger nuclear systems 56 Fe and 
58 Fe, indicating a significant mass number dependence 
of a. We therefore evaluate C a {rj u ) for a ranging from 4 
to 7. 

Figure 3 shows C ' a (Vv) as a function of r\ v for var- 
ious values of a, indicating that the rates of the v- 
process can vary substantially with the adopted values 
for these two key parameters. The production of 7 Li and 
n B is proportional to the reaction retes of 4 He(i/, v'p) 3 H 
and 12 C(i/, i/p) n B, which have similar values of a (see 
above), so that the ejected masses of n B and 7 Li in the 
case of r\ v — 3 would be increased by about 50% in com- 
parison to the yield obtained for r\ v = 0. 

When we allow for nonzero chemical potentials, the 
corresponding range of T v T derived from the GCE con- 
straint for 11 B changes. Consider the relation between 
the neutrino temperatures derived from a given yield ob- 
tained with either nonzero chemical potential T v or with 
zero chemical potential T„o by enforcing A(T„,?7„;i) — 
A(Ti,o, 0; t). The ratio of these two temperatures is given 
as 

^ = cm-^-v. (8) 

For nonzero chemical potentials, T^/T^q is a monoton- 
ically decreasing function of r\ v . In the case of r\ v = 3 
we find T v /T v0 = 0.90 for a = 4 and 0.94 for a = 7. 



This implies that the neutrino temperature satisfying the 
GCE production constraint of n B is reduced to 4.3 MeV 

< T v ^ T (r\ v = 3) < 5.9 MeV, about 6% - 10% smaller 
than the range inferred for r\ v = [Eq. ©]. Likewise, 
T Va and T Pc would be reduced by a comparable fraction. 
In the case of negative r\ v , T v jT V Q increases very weakly, 
e.g., T v /T v0 = 1.015 for j] v = -3. 

In summary, the ejected masses of 11 B and 7 Li in- 
crease with v^ r and v^ r temperature through the en- 
ergy dependence of the cross sections of the ^-process: 
This dependence (oc )is stronger than the depen- 

dence on the total neutrino energy. To reproduce the 
supernova contribution of n B within the framework of 
GCE, neutrino temperature is constrained to 4.8 MeV 

< T v (r) v — 0) < 6.6 MeV. Nonzero neutrino chemical 
potential leads to a larger light element yield. The ejected 
masses of n B and 7 Li would be increased by about 50% 
in the case of rj^ = 3. For a given yield, the required neu- 
trino temperatures are reduced correspondingly, but the 
change is less than 10%. The inferred temperature range 
provides a constraint on theoretical models of neutrino 
transport in supernovae and constrains their 7 Li yields, 
which imposes constraints on contributions from AGB 
stars and novae to Galactic 7 Li. 
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